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Abstract. When a Brownian object is in a nonequilibrium steady state, actual 
force exerted on it is different from one in a thermal equilibrium. In our 
previous paper [Phys. Rev. Lett. 108 (2012), 160601] we discovered a general 
principle which relates the missing force to dissipation rates through a concept 
of momentum deficit due to dissipation (MDD). In this article, we examine 
the principle using various models based on hard disk gases and Brownian 
pistons. Explicit expressions of the forces are obtained analytically and the 
results are compared with molecular dynamics simulations. The good agreement 
demonstrates the validity of MDD. 
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1. Introduction 

Since its inception, understanding the effect of environments on a closed system has 
been a key subject of thermodynamics. In particular when the system size is reduced 
to a mesoscopic scale, fluctuations in the system caused by those of the environments 
play dominant roles in many physical phenomena, such as Brownian motion. The 
state of Brownian objects is typically investigated with the Langevin theory in which 
the environments exert forces on the Brownian objects through deterministic linear 
friction and fluctuating Langevin forces pQ . This approach has been proven to be very 
effective in many applications. Furthermore, the recent development of stochastic 
energetics [2] allows us to investigate the exchange of energy between the Brownian 
objects and the environments within the Langevin theory. 

When the system is in contact with more than one environments which are not 
equilibrium with each other, we expect that energy and momentum flows between 
the system and the environments change in the way that the detailed balance is 
broken. Rigorously speaking, this loss of detailed balance brings the environments 
out of equilibrium at least in the vicinity of the system-environment interfaces. This 
aspect is not reflected in the Langevin description. Consider, for example, cases 
where a Brownian object is simultaneously in contact with two different heat baths 
at different temperatures. What force will be exerted on the Brownian object by the 
baths? A natural extension to the standard Langevin theory is to use the frictions and 
the stochastic forces from each bath assuming that the fluctuation-dissipation relation 
of second kind (i.e. the Einstein relation) holds independently for each bath, the 
condition which is sometimes called local detailed balance. While such a simple linear 
model works well in many cases, there are phenomena that refuse to be understood 
by the linear Langevin approach. 

A most striking example is the adiabatic piston placed between two gases with 
different temperatures [3j [4] . The piston has no internal degree of freedom so that no 
heat flows between the baths through it. An interesting question is if the piston moves 
when the two baths have the same pressure. Naively one may think that the piston 
does not move because pressure on the both sides of the piston is the same. It turns 
out that the laws of thermodynamics alone cannot tell whether the piston moves or not 
[5]. Feynman [6] pointed out that the fluctuations of the piston's velocity should be 
taken into account. However, the Langevin approach with the linear frictions falsely 
predicts zero mean velocity. 

A similar difficulty also appears in some models of Brownian motors working 
between two baths [7 . In these models, the body of Brownian objects, e. g. a 
triangular body, is not symmetric with respect to space inversion. This kind of 
asymmetry is not fully realized in the linear Langevin equation since the linear 
friction constant or tensor is non-polar. In the case of the adiabatic piston, the 
Brownian object itself is symmetric but the environments are not. In either case, 
the linear Langevin theory cannot take into account asymmetric interactions between 
the Brownian objects and the environments. 

A common solution to these problems has been to resort to full and general 
microscopic descriptions, such as molecular dynamic (MD) simulation or master- 
Boltzmann equations under pertinent perturbative approximations. These methods 
are effective in predicting the outcome. For the adiabatic piston, the MD simulation 
[1] and the perturbative mastcr-Boltzmann equation [3J H] give quite consistent 
results showing that the piston moves towards the hotter reservoir. A more recent 
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investigation based on the nonlinear Langevin equation also confirms it [5] . While we 
now know that the piston moves, we still don't fully understand the physics behind it. 

Recently, we investigated the force exerted by gas particles on a Brownian object 
in a non-equilibrium steady state (NESS) and discovered a rather general principle [5] . 
When there is energy dissipation, net momentum flux at the surface of a Brownian 
object is reduced from the flux without dissipation, which we shall call momentum 
deficiency due to dissipation (MDD). As a consequence the force on the Brownian 
object decreases from the equilibrium force by 

7 0) 

Fmdd = -c^ (1) 

Vth. 

where J^. s is energy dissipation per unit time, and the thermal velocity of the gas 
particles of mass m at temperature T is defined by Uth = \J k^T/m. The Boltzmann 
constant is denoted with k^. The positive prefactor c depends on the detail of the 
system but usually at the order of 1. 

With this new principle, we are able to explain all of the above mentioned 
phenomena without the lengthy calculation [S] . It is this MDD that is what is missing 
in the linear Langevin theory. The stochastic energetics [2] tells us that the linear 
Langevin theory is sufficient to obtain the dissipation rate. Therefore, one can evaluate 
the missing force (JXJ) within the Langevin description. Furthermore, we note that this 
fundamental principle is applicable to systems beyond the regular Brownian objects, 
such as inelastic pistons and granular Brownian ratchets (TOl E] • 

In this paper, we will demonstrate the validity of the new principle ([IJ using hard 
disk systems including MD simulation. In the next section, we heuristically explain the 
idea of momentum deficit due to dissipation. Then, we derive explicit expressions of 
-Fmdd for hard disk systems. The results are compared with MD simulation of shared 
Brownian pistons. We also derive non-equilibrium forces on an inelastic Brownian 
piston. Despite that the origin of dissipation is quite different from the previous 
model, we arrive at the same expression, ([I]) and MD simulation confirms it. 

2. Momentum deficit due to dissipation: A heuristic argument 

In this section, we briefly summarize the key concept developed in Ref. [9]. First, we 
consider a simple equilibrium system illustrated in Fig. [lja). The system consists of 
a two-dimensional cylinder filled with a gas and a piston of mass M with surface size 
L, pressed with a constant external force F cxt - The piston is a Brownian object and 
its velocity fluctuates due to the collision with the gas particles. However, when the 
piston is in thermal equilibrium, its mean velocity is zero. Hence, the sum of the pre- 
collisional and post-collisional momentum flows, ^mv^ + ^2 (~ "m-Vout), is balanced 
by the external force. Here the summation is taken over all collisions during a unit 
time. The detailed balance tells us that at equilibrium the two momentum flows must 
be equal and therefore, on average, 2j2 mv 'm = \F cx t\ = pL, where p is the pressure 
of the gas. Note that unlike a simple kinetic theory used in elementary textbooks 
the individual collisions can transfer energy and momentum between the gas and the 
piston at the microscopic time scale since the piston is a Brownian object. It is the 
detailed balance that makes the two momentum flows identical on average. 

Now we turn to a non-equilibrium case shown in Fig. |T|(b) where energy flows 
from the gas through the piston into another environment. We assume that the piston 
is in a NESS so that its mean velocity is zero. The magnitude of the external force is 
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(a) Equilibrium 




(b) Non-equilibrium 




Figure 1. Simple models illustrating the key point of MDD. The upper panel 
(a) shows an equilibrium case where a gas is confined by a piston pressed with a 
constant force F CX f In this case, the momentum flow of the outgoing particles, 
y"](— mn on t) equals the momentum flow of the incoming particles ^2mvi n . The 
sum of the two flows is the pressure of the gas which balances with the external 
force. The lower panel (b) illustrates a non-equilibrium case where the gas loses 
energy through the piston to external environments. In this case, the momentum 
flow of the outgoing particles is smaller than that of the incoming particles, 
resulting in the reduction of the force exerted on the piston by the gas. 



the same as that in the equilibrium case, i. e., _F ex t = pL. In non-equilibrium cases, 
this relation does not necessarily indicate mechanical equilibrium. Even for a non- 
fluctuating macroscopic object, the actual force exerted by the gas deviates from pL 
as observed in a radiometeJJ Our question is what is the actual force on the piston 
when it is a Brownian object. We expect that the post-collisional speed is smaller than 
that in the equilibrium case. Hence, the net momentum flow must be reduced. This 
is the momentum deficit due to dissipation and the force induced by MDD, -Fmdd , is 
defined by 

Fmdb = rnv in + y^-rmw) - pL , (2) 

which vanishes in the absence of dissipation. 

The magnitude of Fmdd can be estimated from the energy balance 

^ ~2 V in ~ "^" W out — -'diss • 00 

As a rough estimate, we replace fluctuating quantities with typical values; v- ln ~ 
lu co \ v t h and v i n ~ w coi ^th w here w co i is the number of collisions per unit time. 
Similarly for the outgoing particles, we introduce a typical velocity w ou t- Furthermore, 



| Maxwell struggled to explain it until near the end of his life. See Ref. 1121 . 
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(a) Shared Brownian piston 






• •* 

*. . ' . *. 








• • . m 

•• • 

• T*Pt . • * . 
• • • 

• • • • 


• ••••• »• 



(b) Inelastic Brownian piston 
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Figure 2. Models with two different types of dissipation, (a) Shared Brownian 
piston: The piston is in contact with the second gas. When two gases have 
different temperatures, heat flows from the upper gas to the lower gas through 
the fluctuation of the piston, (b) Inelastic Brownian piston: The collision between 
the gas particles and the piston is inelastic. The energy dissipates into the internal 
degrees of freedome in the piston and the gas particles. 



the dissipation is assumed to be so weak that i>i n — u ou t ~ 2w t h- Then, the balance of 
momentum and energy are expressed in simpler forms: 

Fmbd ~ -Wool m(v t h + v out ) (4) 

^dii ~ w ™i m(v t h + Vout) v t h ■ (5) 

Eliminating the unknown quantity V ou t, we obtain Eq. ([1]) except for the prefactor 
c which is omitted in the above phenomenological argument since it depends on the 
system configuration. 

Despite the drastically simple derivation, Eq. (JTJ agrees with the result of lengthy 
calculations except for the prefactor. Applying this principle we are able to explain 
the adiabatic piston and other phenomena both conceptually and quantitatively [5]. 

3. Momentum deficiency due to dissipation in hard disk systems 

3.1. Basic model 

We consider again the NESS case shown in Fig. Hfb). In order to find explicit 
expressions we assume that the gas consists of hard disks which elastically collide 
with the piston. The temperature of the gas, T, is assumed to be constant and 
the velocity of the gas particles satisfies the Maxwellian velocity distribution^ We 

§ In general, this assumption does not hold under nonequilibrium conditions. Following the previous 
models 4, 7 10 IT], we assume that the incoming particles leave a thermostated region and directly 
hit the Brownian object. The outgoing particles, lower in kinetic energy on average, travel back to 
the thermostated area before colliding the incoming particles, and thus the incoming particles follow 
the Maxwellian distribution. However, this situation is possible only when the distance between the 
Brownian object and the thermostated region is not larger than the mean free path. 
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assume only binary hard collisions to take place. The piston surface is smooth so that 
the velocity component parallel to the piston remains the same upon the collision. 
Hereafter, we consider only the velocity component perpendicular to the piston. 
Momentum and energy conserve at each collision even when the Brownian object 
is simultaneously in contact with other baths or external agents since the hard disk 
collision is instantaneous. For i-th collision, the gas particle and the piston have the 
pre-collisional velocity Vi and Vi, respectively and the corresponding post-collisional 
velocities Ui and Ui are determined by the momentum and energy conservation laws: 

mil, - mvi = MVi - MU, (6) 
m o m n M 9 M 

~2 1 ~ ~2 1 = ~J * ~ ~2 * ' (7) 
Between the successive collisions, the piston interacts with another environment or 
an external agent and its momentum and energy change by APt and AEi, respectively. 
The change in the piston velocity is determined by another set of momentum and 
energy balance equations: 

MVi+i - Am = AP, (8) 
M 2 M , 

Y 1+1 ~ T 

Summing up Eqs. © and ([5]) over all n collisions during a unit time, we find the net 
momentum balance: 



i ^V& 1 -—Uf = AE i . (9) 



Wcol (m(v) col - m(u) col ) + F cxt = , (10) 

where F cx t = Yl an d we have used the NESS condition MU n — MV\. The mean 
value is defined by J^Vi — w co i(v) co \. Note that (■ ■ -) co i indicates average over all 
collisions during a unit time (see Appendix) and it is not the same as a regular thermal 
average over the Maxwell distribution. The force due to MDD is now expressed as 

-Fmdd = Fcxt -pL= -w co i m {{v) co \ + {u) co \) . (11) 
Similarly, adding up Eq. (J7J along with Eq. ^ leads to net energy balance 

4i = ^(^(u'U-^(v 2 ) col ) , (12) 

where J^ ss = £ AE l . We have used the steady state condition MUl/2 = MV?/2. 

For incoming particles, we find (v) co \ — ^/7r/2u t h and (v 2 ) co \ / '(v)l ol — A/it (see 
Appendix). We do not have exact statistics of the outgoing particles. However, when 
the dissipation is weak, we can assume that (u) co \ ~ — {v)co\ and (w 2 )coi/(w)coi ~ 
( u2 )coi/(u)l ol - Using this approximation, Eq. (|T2^) is reduced to 

^diis = «thWcol m {(v) co i + (u)coi) ■ (13) 
Comparing Eqs. (fTTj) and (|T3|) . we obtain the formula ([lj with a prefactor c = W7r/8. 



3.2. Shared Brownian pistons 

Now, we introduce a more concrete NESS model shown in Fig. [Ha) so that we can 
evaluate the dissipation rate. The upper cylinder is the same as the basic model [Fig. 
Hfb)] and the gas in it has temperature T\. A NESS condition is generated by linking 
the piston to another piston in the second cylinder filled with another gas at a different 
temperature T^. These two pistons are rigidly connected and move together. Unlike 
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Figure 3. Molecular dynamics simulation of momentum deficit due to dissipation 
using the model illustrated in Fig. EI a). The upper panel shows the position of the 
piston and the lower panel the temperature of two gasses. Initially,the two gases 
have the same temperature T\ = T2 = 1.0. When the temperature of the lower 
gas is reduced to T2 = 0.5, the upper gas is compressed due to Fmdd- The values 
of the system parameters are L = 300, M/m = 20, Xq = 1000, pi = P2 = 0.00333 
where the diameter of the gas particle is a unit of the length. The mean free path 
in this configuration is about 150. The data are averaged over 300 realizations. 



the upper one, the lower cylinder is periodic so that the particle density does not 
change as the piston moves. Accordingly, when the piston moves, the pressure of the 
upper gas, p\ changes while the pressure P2 of the lower gas remains constant. 

It is known that when T± > T 2 , heat flows from the upper gas to the lower gas 
through the fluctuations of the Brownian object [13] ■ The heat dissipation through 
a shared piston is understood at the level of standard Langevin theory [14]. Using 
the friction coefficients for the upper and lower pistons, 71 = 
% j -k p 2 L\/ k B T 2l we find the dissipation rate as 

7r k B Ti - k B T 2 



'8/tt p\L\fkoT{ and 



72 



J. 



(e) 



s M( 7 r i + 72 - 1 )- 

Substituting this dissipation rate and the prefactor c 
an explicit expression of the force due to MDD: 

2pip 2 L m 



MDD 



(fc B Ti - k B T 2 ) 



(14) 

to Eq. (fTJ , we obtain 
(15) 



pi + 2p 2 M 

We have checked the above results using hard disk molecular dynamics simulation. 
The detailed simulation method will be written somewhere else. Initially the system 
is at a thermal equilibrium with T\ = T% = 1.0. The mean position of the piston X 
remains constant since F cxt + PiL — 0. Figure [3] shows that when the temperature 
of the lower gas is reduced to T 2 = 0.5, the upper gas is compressed despite the 
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Piston mass [M/m] 

Figure 4. Molecular dynamics simulation of momentum deficit due to dissipation 
using the model illustrated in Fig. [2j a )- The black solid circle shows the force 
measured in the simulation through Eq. I|16|l . The red solid square plots the 
force estimated from Eq. JTJ using the observed heat flow. The full theoretical 
prediction Eq. I|15| l. is plotted with a solid line. The upper curves shows the 
case where the temperature of the second gas (T2 = 1.5) is higher than the upper 
gas whereas the lower curves shows the opposite case where the lower gas has 
T2 = 0.5. See Fig. [3]for the values of the system parameters. 



temperature is kept at T\ = 1.0. The displacement of the piston indicates that the 
force exerted on the piston by the gas is not the pressure times the surface area. 
Similarly, when T2 is raised above Ti, the upper gas expands. 

When the system reaches a NESS, the piston is settled at a new position and a 
new pressure p\ is established. Assuming that the gas obeys the ideal gas law and the 
displacement AX is much smaller than Xo, the missing force is estimated by 

AX 

Fmdd = {pi - p[)L = PiL -— . (16) 

In Fig. @] we plot Fmdd obtained in three different ways; Eq. (| 10[) with the 
measured displacement of the piston, Eq. ([1]) using the dissipation rate measured in 
the MD simulation, and the full theoretical result (fl5l) . All three estimations agree 
very well, implying the validity of Eq. (UJ). 



4. MDD in granular systems: Inelastic piston 

In order to demonstrate the generality of Eq. ([1]), we consider a different type of 
dissipation. In the model illustrated in Fig. [U[b) energy dissipates into the internal 
degrees of freedom of the piston and gas particles though inelastic collisions between 
them. Using a standard collision rule used for granular systems, post-collisional 
velocities are related to the pre-collisional velocities as it, — 11% = —e(vi — Vi) where e 
is a coefficient of restitution (0 < 1 — e 1). The momentum always conserves and 
thus Eq. © is still valid. The energy balance for this model is 

m 9 m ? M ~ M rr n 1 — e 2 vnM , 2 

2 vl ' = ~2 1 - t 1 + —^nt^ - ^ (17) 
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Figure 5. Fmdd on the inelastic piston shown in Fig. [2fb) is plotted as a 
function of the piston mass (upper panel) and of the restitution coefficient (lower 
panel). The force was evaluated in three different ways; the black circle indicates 
the force measured from the displacement of the piston in the MD simulation 
using Eq. Ill6jl . The red square shows the formula {TJ using the dissipation rate 
observed in the MD simulation. The solid line plots the theoretical value, the 
sum of Eqs. (1_'21 and l|23|l . In the upper panel, e = 0.96 and in the lower panel, 
M/m = 20 are used. See Fig. [3] for other parameter values. 



Unlike the previous case, there is no dissipation between collisions, hence V^+i = E/j. 
Summing up Eq. (|17j) for all collisions during a unit time, we obtained Eq. (|12|) again 
with a different dissipation rate: 



I I 1 III M , 

JSL == 2 m + M Ucol{{V " V) )co1 (18) 

7(c) 



Although the actual expression of JjJL 1S different, the momentum and energy 
conservation laws are universal, and hence the general expression (fTJ) is also valid 
for this model. 

Now, we evaluate the dissipation rate (TT51) . Assuming that the piston obeys the 
Maxwellian distribution with a kinetic temperature Tki n , the mean value in Eq. (|18[) 
is given by 

i, m 2\ 2k B T 2fc B r k i n , , 

((« - V) coi = + — — — , (19) 

m M 
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as shown in Appendix. The dissipation rate is divided into two parts, one proportional 
to the mean kinetic energy of the gas particles (1st term) and the other to the mean 
kinetic energy of the piston (2nd term). We shall call the former house keeping 
dissipation (Jdiss,hk) and the latter excess dissipation (Jdiss,cx) as coined in Ref.[15 . 
Using the mean values used in the previous section and w co i = pLvth/v / 27r (see 
Appendix), we obtain 

Jd iss.hk 

ISS.CX 

where we assumed m/M < 1, (1 - e 2 ) « 2(1 - e), and T kin w T X (1 + e)/2 [16]. 
Through our basic principle (QJ and the prefactor c = these dissipation rates 

lead to Fmdd = i*MDD,hk + i*MDD,ex where 

i J1 MDD,hk = -^(l-e)pL (22) 

Tfi 

fjvlDD.ex = -77(1 - e )P L ■ (23) 

Figure [5] shows the result of MD simulation. We again plot -Fmdd evaluated in 
three different ways; Eq. (fl6|) . Eq. (fT|) using the measured dissipation rate, and 
full theoretical value with Eqs. (|2"2l and (|2"3")l . All three estimations agree well. 

As discussed in Ref. [9] , these forces explain the driving of inelastic pistons [10] 
and granular ratchets [TT]. For the granular pistons, the house-keeping force (|2"2"j) 
is dominant and agrees with the perturbative results |10) . On the other hand, the 
driving force of the granular ratchets is the excess-dissipation force (|23[) since the net 
house-keeping force vanishes in this model. 

5. Conclusions 

We examined the momentum deficit due to dissipation using various hard disk models. 
The explicit expressions of the force due to MDD are obtained for two different models; 
the shared Brownian piston and the inelastic Brownian piston. Despite the different 
dissipation processes, the general principle ([T]) is valid for both cases. Molecular 
dynamics simulations agreed well with the theoretical predictions for both cases. 
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Appendix A. Collision statistics 

Here, we briefly explain the calculation of average over collision events. We assume 
that the incoming gas particles obey the Maxwell's velocity distribution f g (v) = 
y/m/2irk B Te- mv / 2ksT with temperature T. Under non-equilibrium conditions, the 



(1 - e)\j —vthpL 



(l-e) 



7 k B T 
M 2 



(20) 
(21) 
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velocity distribution of the piston is not necessarily Maxwellian. However, when the 
piston is under a NESS not far from equilibrium, the piston approximately follows 
the Maxwell's velocity distribution f p (V) = v / M/2TTk B T kin e~ Mv2 '/ 2k ^ T ^ but with a 
kinetic temperature T k { n which is different from the temperature of the gas. Under 
these assumptions, the velocity distribution of gas particles colliding with the piston 
moving at velocity V is given by 

$(„; V ) = Lp{v ~ V) f g (v)f p (V)e(v - V) (A.l) 

Wcol 

where 0(-) is the Heaviside step function and the normalization constant lj co \ is the 
total number of collisions per unit time defined by 

/>oo />oo 

Wool = Lp dV dv(v - V)f g (v)f p (V) 



= L JMJ 1 + ^,i^ (A.2) 

where m <§; M is assumed. 

Using the probability distribution (|A.1[) we obtain the 1st and 2nd moments: 

/>oo />oo 

(»)col = / dV dvv$(v;V) 



irk B T mT kin \ 2 hr 

^ T7^r~ ~\h V th (A.3) 



2vl (A.4) 



2m \ MT r 

/oo poo 
dV dvv 2 ${v;V) 
-oo J —oo 

_ keT mT kin + 2MT 
m mT k i n + MT 
Similarly, the second moment of the relative velocity is computed as 

/oo y* oo 

dV dv(v-V) 2 <S>(v;V) 
-oo J —oo 

= 2k B T 2fc B T kin 
m M 
which is exact under the present assumption. 



(A.5) 
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